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I. INTRODUCTION 



The search for the Higgs boson will be a top priority of the 
CERN Large Hadron Collider (LHC). The LHC's discovery 
potential for the Standard Model Higgs boson fully covers the 
mass range from the experimental lower bound established by 
the LEP experiments (M H > 114 GeV) up to M H ps 1 TeV, be- 
yond which the concept of the Higgs boson as an elemen- 
tary particle becomes questionable. In addition, the Fermilab 
Tevatron could find evidence for or even discover the Higgs 
boson if Mh ^ 180 GeV and if sufficient luminosity yl can 
be collected. 

The theoretical description of the signal processes for Stan- 
dard Model Higgs boson production is under good control. 
For a review, see Ref. |2]. The dominant production mode 
is gluon fusion, for which the next-to-next-to-leading order 
(NNLO) corrections are now available |3, 4] and have recently 
been reconfirmed by Ref. 1 5 ] . The radiative corrections for the 
weak boson fusion channel 1 6 ] and associated production with 
a weak gauge boson @] have been known for several years, 
rendering the theoretical uncertainty in these processes very 
small. Recently, next-to-leading order (NLO) corrections have 
also been evaluated for Higgs boson production in association 
with top quarks IMI^ Ilfllllll . resulting in a drastic reduction 
of the scale uncertainty. 

These results can be used for supersymmetric Higgs boson 
production as well. However, because of the enriched par- 
ticle spectrum in supersymmetric extensions of the Standard 
Model, they provide only a part of the full production rate in 



general. Additional contributions arise through intermediate 
supersymmetric partners fl2Tl and modified couplings of the 
Standard Model particles. In order to avoid unnecessary gen- 
eralizations, we will focus on the Minimal Supersymmetric 
Standard Model (MSSM) for the rest of this paper (see, e.g., 
Ref. IL3I1 for an outline of the MSSM). The extent to which 
our results can be transferred to other models should be clear 
from this discussion. 

The MSSM contains two Higgs doublets, one giving mass to 
up-type quarks and the other to down-type quarks. The associ- 
ated vacuum expectation values are labeled v„ and vy, respec- 
tively, and they determine the MSSM parameter tan j3 = v u /Vd- 
After spontaneous symmetry breaking, there are five physical 
Higgs bosons, whose mass eigenstates are denoted by h ("light 
scalar"), H ("heavy scalar"), H ± ("charged scalars"), and A 
("pseudoscalar"). One interesting consequence of this more 
complicated Higgs sector is that, compared to the Standard 
Model, the bottom quark Yukawa coupling can be enhanced 
with respect to the top quark Yukawa coupling. In the Stan- 
dard Model, the ratio of the ttH and bbH couplings is given at 
the tree-level by A r SM /A^ M = m t /mj, w 35. In contrast, in the 
MSSM, it depends on the value of tan /3 . At leading order, 
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where a is the mixing angle between the weak and the mass 
eigenstates of the neutral scalars. A value of tanj3 as large as 
30 — 40 could be accommodated fairly naturally in the MSSM. 
Such an enhancement would have (at least) two important 
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FIG. 1 : For large tan /3 , the bottom quark contribution to the gluon 
fusion process can be comparable to the top quark contribution. 



consequences. The first is that in the gluon fusion mode it 
is no longer sufficient to consider top quark loops as the only 
mediators between the Higgs boson and the gluons; one must 
also include the effects of bottom quark loops (see Fig.Q. 
Since one cannot justify the use of an effective field theory in 
which the bottom quarks are integrated out, this involves com- 
puting massive multi-loop diagrams. While the massive NLO 
calculation (including massive two-loop virtual diagrams) was 
performed some time ago 0, the NNLO result (requiring up 
to three loops for the virtual correction) is still beyond the lim- 
its of current calculational technology. 

The second consequence is that Higgs boson production in as- 
sociation with bottom quarks can become an important chan- 
nel: pp — > bh(j) (p £ {p,p} and <p 6 {h,H,A} here and in 
what follows). At first sight, the evaluation of the correspond- 
ing cross section is in close analogy to the process pp — > tt<f>. 
But this is only true if the bottom quarks are observed in the 
detector, and are thus restricted to large transverse momenta. 
If at least one of the bottom quarks escapes detection, the pro- 
duction rate must be integrated over all transverse momenta 
of this bottom quark. Since the Higgs boson is much heavier 
than the bottom quark, this integration leads to collinear log- 
arithms, \n{tnb/Mty), which require a more careful analysis 
than in the case of tt production. 

The subject of this paper is the precise prediction of the to- 
tal cross section for Higgs boson production in association 
with bottom quarks, where neither bottom quark need be de- 
tected. This requires integrating over the transverse momenta 
of both final state bottom quarks. Each integration gives rise 
to collinear logarithms of the kind mentioned above. Since the 
bottom quarks may remain undetected, it is more appropriate 
to view our result as a part of the total inclusive Higgs produc- 
tion rate o(pp — > <j> +X). In order to emphasize this point, we 
shall henceforth denote the fully inclusive process mediated 
through bottom-antibottom annihilation as pp—> (bb)^ + X. 

Our calculation is based on the approach of Refs. 1 15, 
where the leading-order (LO) partonic process is taken to be 
bb — > (j). We will refer to this as the variable flavor num- 
ber scheme (VFS) approach in what follows, as opposed to 
the fixed flavor number scheme (FFS) approach, where the 
tree-level process gg — ► bb§ is taken as the lowest-order 
contribution and bottom quarks cannot appear in the ini- 
tial state. The initial state bottom quarks in the VFS ap- 
proach arise (predominantly) from gluon splitting in the pro- 
ton, parametrized in terms of bottom quark parton distri- 
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this way, the large collinear logarithms that arise due to the 
fact that the colliding gluons carry a momentum of the or- 
der of M§ /2 ^S> nif, can be resummed through DGLAP evolu- 
tion. The convolution of these bottom quark densities with 
the partonic cross section leads to the hadronic cross section 
o[pp->(bb)<j>+X]. 

This process has been a subject of interest for some time. 
It is currently known up to NLO in the VFS approach fl5l 
E [13 EH 0. In the FFS approach, the calcula- 
tion is analogous to tt (j) production, which is also known to 
NLO H H [jjj El EH lH H. The case where one bottom 
quark is tagged has been computed at NLO in Ref. |34[; in 
that case, the LO process in the VFS approach is bg — > b(p. 

There has been an ongoing discussion as to the relative merits 
of the VFS and FFS approaches flfil QJ, | H l3fil l3l . es- 
pecially because the results of the two approaches disagree 
by more than an order of magnitude for the scale choice 
= M« = Mq, where jip and /Ir are the factorization and 
the renormalization scales, respectively. Recently, it has been 
argued 1 17], that the proper factorization scale for this process 
should be /Xf rs Ma /4 instead of Mq. Indeed, for this choice, 
the disagreement between the VFS and the FFS approach is 
significantly reduced. The result of our paper demonstrates 
the self-consistency of the VFS approach and confirms the 
proposed factorization scale of Ref. 11711 as the appropriate 
choice. 

Given the considerations above, the motivations for a NNLO 
calculation are manyfold. One is to examine the assertion 
that Uf ~ Mq/4 is the proper choice for this process at 
NLO OS. If the higher order corrections are minimal at 
that scale, this would be a strong argument in favor of the va- 
lidity of the VFS approach to (bb)<j> production. A second mo- 
tivation, as will be discussed in more detail below, is that the 
NNLO terms play an exceptional role in the VFS approach to 
(bb)ij) production due to the fact that they are the first to con- 
sistently include the "parent" process, gg — > bixj), and thereby 
sample the same range of bottom quark transverse momenta 
as the LO FFS approach. A third and perhaps dominant moti- 
vation for the NNLO calculation is to reduce the sensitivity of 
the calculation to the unphysical scale parameters jif and 
thereby removing a significant source of uncertainty from the 
theoretical prediction. 

In this paper we will present results for the process pp — > 
(bb)§ +X at NNLO. As will be shown, they nicely meet all 
expectations concerning their dependence on the renormaliza- 
tion and factorization scales, thus providing a solid prediction 
for the total cross section at the LHC and the Tevatron. The 
inclusive production cross section could have phenomenolog- 
ical implications for the observation of the supersymmetric H 
and A bosons, for example, in the H/A — > decay mode. 

The organization of the paper is as follows: In Sec. [H] we dis- 
cuss the VFS approach to computing the pp—> (bb)(j) +X pro- 
cess and its motivations. In Sec. H!]]we describe the actual 
calculation and in Sec. IIVI we present our numerical results. 
Analytic results for the partonic cross sections are presented 
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in the appendix. 
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FIG. 2: Partonic processes for pp — * bbH. Not shown are diagrams 
that can be obtained by crossing the initial state gluons, or radiating 
the Higgs off an antibottom quark. 



II. THEORETICAL DESCRIPTION OF THE 
PRODUCTION RATE 



In the FFS approach, the LO partonic process for the produc- 
tion of a Higgs boson in association with a bottom quark pair 
is of order a 2 . A few typical diagrams are shown in Fig.|2] Be- 
cause of the large mass difference between the bottom quark 
and the Higgs boson, the total cross section contains large log- 
arithms of the form 



(2) 



where jj,^ is of the order of M§. More precisely, every on- 
shell gluon that splits into a bb pair with an on-shell bottom 
quark generates one power of that logarithm. Thus, Figs.|2ja) 
and|2jb) generate two and one power of 4, respectively, while 
Figs.|2jc) and |2jd) do not generate any 4 terms. Further- 
more, each higher order in perturbation theory brings in an- 
other power of If, due to the radiation of gluons from bottom 
quarks. 

Because ly ~ \n(mj y /M^) is rather large, aJb is not a good ex- 
pansion parameter. It would be better to re-organize the per- 
turbative series such that terms like (a s lb) n are resummed to 
all orders in n. This resummation can be achieved by introduc- 
ing bottom quark parton distribution functions which contain 
all the collinear terms arising from the spl i tting of gluo ns into 
bb pairs JIB. 113 EH 13 E3, H3, 13 M Eli. This 
constitutes the motivation for using the VFS approach. 

Convolving the tree-level process of Fig.[3la) with these bot- 
tom quark distributions will resum the leading logarithms of 
the form (oc s lb) 2 ■ (a s lb) n , n > 0. In order to retain subleading 



logarithms, one has to compute higher orders. For example, 
including the NLO contributions with all the relevant subpro- 
cesses {bb — > h, bb — ► hg, gb — > hb, and gb — ► hb), resums the 
terms of order a 2 lb (of s /&)", n > 0. In order to retain all pow- 
ers of lb at order a 2 £„(o£A)", it is necessary to evaluate the 
cross section up to NNLO. 

Let us briefly review the idea of the VFS in its simplest form. 
Assume ng =rif—\ massless quark flavors and one massive 
quark flavor of mass ;«/,. First, one defines parton densi- 
ties fj"^ (x, Q 2 ) for the gluon (i = g) and the massless flavors 
(i = 1, . . . , n£) in the standard way, obeying DGLAP evolution 
with rii active flavors. The heavy quark density fn" e \x,Q 2 ) 
is assumed to vanish. Partonic processes involving the heavy 
quark should be evaluated by keeping the heavy quark mass. 
This is called the n<>-flavor scheme. 

At a certain scale \xi, one relates the n#- to the n^-flavor 
scheme by defining initial conditions for new parton densities 
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The Qj are determined by the requirement that physical quan- 
tities are the same (up to higher orders in a s ) in both the rip- 
and the n/-flavor scheme. (This requirement may be imple- 
mented asymptotically or using mass dependent terms l20tl2ll 
1241. l25ll .1 Above the matching scale, the DGLAP evolution of 



thcfj" f \x,Q 2 ) 
flavors. 



(2 = g, 1 , . . . ,nf) is performed with n f active 



In general, one assumes the n^-flavor scheme at scales Q 2 < 
m? and switches to the n ^-flavor scheme at larger values of Q 2 . 
It is also convenient to choose the matching scale in Eq. © 
as jl^ = m 2 h , which avoids the occurrence of logarithms of 
m h /n h . 

For our purposes, mi 

IQ 2 ~ m 2 /M 2 < 0.003, so threshold ef- 
fects from the matching prescription should be minimal. For 
the same reason, it is justified to neglect the bottom quark 
mass in the partonic process (apart from Yukawa couplings, of 
course). Indeed, masses must be neglected for initial state bot- 
tom quarks in order to avoid violation of the Bloch-Nordsieck 
theorem (and the consequent failure to fully cancel infrared 
divergences) at NNLO and beyond 1 39]. 

In order to make the following discussion more transparent, 
let us write the fully inclusive (bb)<j> production rate in the 
VFS approach schematically as follows: 

o(pp^{bb)$+X) = 

£ (aJi,)" \ a, 2 [cm if, + c n\ h + cm] 



The sum over n is implicit in the parton densities. A LO cal- 
culation determines the coefficients c„q, while NLO adds the 
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FIG. 3: Lowest order diagrams contributing to (a) bb — » H, (b) bb — » Hg, and (c) fcg — > /ffc. At NNLO, these diagrams receive corrections up 
to two loops in case (a), and one loop in case (b) and (c). 



coefficients c n \. Note that the subprocess bg — ► bp does not 
fully determine the coefficients c n \\ in order to obtain the cor- 
rect resummation at a"l'^ 1 (n > 2), one has to include the 
real and virtual corrections to the bb — > p subprocess as well. 
In the same way, the sum of all subprocesses that contribute 
at second order determines c„2 (n > 0), and thus all terms as- 
sociated with the order aj (a s lb) n . The NNLO result is thus 
the first to include all terms of order 0% (as well as higher 
order terms resummed in the parton distribution functions). 
Higher orders in perturbation theory correspond to the coeffi- 
cients £/„(.; their 4 terms are — formally — completely con- 
tained in the parton densities. This illustrates once more the 
exceptional role of the NNLO corrections in this approach. 



The leading order terms were evaluated by Eichten et al. |2S 
The leading bg and gg initiated processes (Figs.|3jc) and|2ja)) 
were subsequently added by Dicus and Willenbrock 1 15|. Ten 
years later, Dicus et al. 11611 (see also Ref. 13511 ') computed the 
full NLO contribution to bb — ► p (and related subprocesses), 
which leads to the single logarithmically suppressed term c n \ 
for all n > 0. Setting the renormalization and the factoriza- 
tion scale equal to {jiR = fif = Mq), they found that the 
ff(<Xs) -corrections to the bb — ► p subprocess and the contri- 
bution from the tree level bg — > pb subprocess are quite large, 
but of opposite sign. This leads to large cancellations which 
are particularly drastic at the LHC. They also observed that the 
contribution from bg — ► (j)b becomes especially large at Higgs 
boson masses below » 150 GeV, meaning that logarithmically 
suppressed terms become important in this region. 

Recently, Maltoni et al. lll7ll revisited the NLO calculation in 
the light of Ref. 1 38], which gives an argument for the proper 
choice of the factorization scale when using the bottom quark 
density approach. Following that argument, they determined 
the factorization scale for the (bb)p process to be [if w M<j, /4. 
With this choice, both the NLO corrections from the gb and the 
bb initiated process turn out to be very well behaved. 

As we will show in Sect. lIVI the behavior of the NNLO correc- 
tions confirms this choice of scale at NLO, in the sense that the 
perturbation theory up to NNLO is very well behaved for this 
choice. This supports the method of Refs. ITvl 13^1 (see also 
Ref. 14011 ') for determining the factorization scale in the VFS 
approach at lower orders. For the process under consideration, 
it turns out that the dependence on the unphysical scales of the 



NNLO result is so weak that the discussion on the proper scale 
choice becomes irrelevant. The overall conclusion is that the 
prediction for Higgs boson production in bottom quark fusion 
is now under very good control. 



III. OUTLINE OF THE CALCULATION 



As discussed before, we will neglect the bottom quark mass 
everywhere except in the Yukawa couplings. The calculation 
is thus completely analogous to, say, Drell-Yan production of 
virtual photons 1 31 14111 : One evaluates virtual and real correc- 
tions to Higgs production in bb, gb, gg, bb, qb and qq scatter- 
ing (and the charge conjugated processes) and then performs 
ultraviolet renormalization and mass factorization. 

The subprocesses to be evaluated at the partonic level are 
given as follows (q S {u,d,c,s}): 

• up to two loops: bb — > p [Fig.|5fa)] 

• up to one loop: bb — > (j)g, gb — » (j)b [Figs.|3jb),|3jc)] 

• at tree level: bb — > <pgg, bb — > pqq, 

bb — > pbb, gb — > (pgb, 

bb -► pbb, bq -> <j)bq [Figs.0a)-Eff)] 

gg^pbb [Figs.aa)-Etc)], 

qq-^pbb [Fig.Hd)] 

We compute the two-loop virtual terms by employing the 
method of Refs. i42l l43ll . which maps them onto three-loop 
two-point functions. In this way, they can be reduced to a 
single master integral, using the reduction formulas given 
in Refs. |44j, 14511 . The master integral has been computed 
in Ref. |46]. The pole parts of this "bbp form factor" can 
be compared to the general formula of Refs. 1 47, 48] which 
provides a welcome check. For the generation of the diagrams 
we use QGRAF 1 49] as embedded in the automated system 
GEFICOM iSPlUlll Vsee also Ref. jH). 

The one-loop single emission processes are obtained by com- 
puting analytically the full one-loop amplitudes, which are 
then interfered with the amplitudes of the tree-level processes. 
The two-particle phase space integrals are also computed an- 
alytically. 
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FIG. 4: Diagrams contributing at NNLO. Note that the Higgs boson can couple to the fe-quarks at any point; only representative diagrams are 
shown. 



For the tree-level double emission processes, we express the 
matrix elements and phase space measures in terms of the 
variable x = / S, where S is the center of mass energy. Then 
we expand the integrands in terms of (1 — x) \3\. This leaves 
us with only one nontrivial phase space integral, independent 
of the order of the expansion. The regular integrands and the 
finite integration region ensure the validity of the interchange 
of integration and expansion. Keeping of the order of ten 
terms in the expansion in (1 — x) leads to a hadronic result 
that is already phenomenologically equivalent to the analytic 
result. By evaluating the expansion up to sufficiently high or- 
ders, however, one can invert the series |53j] by mapping the 
expansion onto a set of basis functions. The latter can be de- 
duced from the known NNLO Drell-Yan result ma 

All algebraic manipulations are performed with the help of the 
program FORM 1541 . 

For a consistent treatment of the NNLO process, it is not suf- 
ficient to evaluate only the partonic cross section at NNLO. 
Another ingredient is the proper parton densities, obeying 
NNLO Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
evolution. At present, only approximate evolution kernels 
are known, derived from moments of the structure func- 
tions 1 551 1561 15711 . On this basis, approximate NNLO parton 
distribution sets have been evaluated 1 58 ] . We use this set in 
all of our numerical analyses below. Once parton distribu- 
tions that use exact NNLO evolution become available, it is a 
straightforward task to update the analysis using the partonic 
results presented in AppendixlAl 

Let us now turn to the underlying interaction and the renor- 
malization of the partonic results. We ignore the bottom quark 
mass and the electroweak interactions, so for our purposes, the 
Lagrangian is: 

= ~F^F a ^+"£qmq + bWb-^b, ( 5 ) 

<7 



where F£ v is the gluon field strength tensor, is the QCD co- 
variant derivative, and the sum runs over the quarks u,d,s,c. 
Xjf is a bare bottom Yukawa coupling constant. In the modi- 
fied minimal subtraction (MS) scheme, the scalar coupling is 
renormalized such that 1 

X b B = X b Z m {a s ) , [<l>=h,H} 

Z„,(a,) = l-^- 
% e 

/a s \ 2 T 1 /15 n f \ 1/101 5 V 

+ ^(a s 3 ), 

(6) 

where e = (4 — D)/2 and D is the number of space-time di- 
mensions in which we evaluate the loop (and phase-space) in- 
tegrals. Z,„(a s ) is identical to the quark mass renormalization 
constant of QCD 1 59, 60]. Here and in the following, we use 

the short hand notations \ b = X^diR) and a s = aj"^(^ fl ) 
for the MS -renormalized Yukawa and strong coupling con- 
stants, respectively. jiR is the renormalization scale, and iif 
is the number of "active" quark flavors. We will set iif = 5 in 
our numerical analyses. 

There are (at least) two methods of obtaining the result for 
pseudoscalar production. The first is to replace the Yukawa 
interaction term in Eq. (|5jl with a pseudoscalar interaction, 

kj?tyb—>ikj?b$Kb, (7) 
and proceed by direct calculation. 



We refrain from quoting terms proportional to |fi and ln47T that drop out of 
MS -renormalized quantities. 
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The second method is to exploit the chiral symmetry of the 
bottom quarks in Eq. which implies that we are free to 
perform independent left-handed and right-handed phase ro- 
tations of the bottom quarks. If we perform the rotation 



bR ib' R 



L i 



(8) 



the Lagrangian becomes 



]_ F a pafiv 

4^ r 



Viipb' 



agb'tYsb', 



(9) 



and we find the same interaction Lagrangian as in Eq. Q. 
This implies that the cross section for pseudoscalar Higgs bo- 
son production, written in terms of the Yukawa coupling Xb, 
is identical to the cross section for scalar Higgs boson produc- 
tion to all orders in <X S . 

Following the prescription of Larin l6lll 2 for the treatment of 
75 in dimensional regularization, we have performed the direct 
calculation through NNLO and find that this is indeed the case. 

Even in the direct calculation, one can see that this identity 
will hold to all orders in a s with the following argument. If 
we square the amplitude before computing loop integrals, all 
fermion lines are closed loops. The fact that we set the bot- 
tom quark mass to zero means that both Higgs boson vertices 
(in both the scalar and pseudoscalar cases) must appear on the 
same fermion line. If only one Higgs vertex were to appear 
on a fermion line, there would be an odd number of y ma- 
trices in the fermion trace which would therefore vanish. In 
the pseudoscalar case, this means that nonvanishing fermion 
traces must contain either zero or two 75 matrices. The pre- 
scription of Larin |61] allows one to assume anticommutativ- 
ity of 75 and identify 7^ = 1 when two 75 -matrices are on the 
same fermion line. Thus, the 75-matrices can be eliminated 
and we see that the calculation for pseudoscalar Higgs boson 
production is identical, diagram by diagram of the squared 
amplitude, to that for scalar Higgs boson production, apart 
from the different Yukawa couplings. 

For the sake of completeness, let us remark that the Stan- 
dard Model value for the coupling constant is given by Xb = 
\flmi, /v, where v « 246 GeV is the vacuum expectation value 
for the Higgs boson field, and m& is the running MS mass of 
the bottom quark, my{iin), evaluated at the renormalization 
scale )1r. In the MSSM we have 



V2- 



nib sin a 
v cos/3 

mb cos a 



v cos/3 ' 
V2^ tan/3, 



=h, 

$=H, (10) 
</> =A. 



The renormalized partonic results have a dependence on the 
unphysical scales /Xp and /Xr , both explicitly in terms of log- 
arithms, and implicitly through the parameters Ok(/i#) and 
A/>(£tfl)- The variation of a s and Xb with is governed by 
the renormalization group equations (RGEs) 



lil—^a s = J3 {a s )a s , fij—jXb = f{a s ) X, 
where 

/3 (a,) = -a s Po - a% - a% + G{a A s ) , 
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Here, = £(n) is Riemann's ^-function (£3 w 1.20206). In 
order to evaluate a s (fJ,R) from the initial value 3 a s (Mz), /3 (a s ) 
is expanded up to af , with I = 1 at LO, I = 2 at NLO, and 
i = 3 at NNLO. The resulting differential equation of Eq. (II It 
is solved numerically. 

In order to evaluate Afc(jU) from its initial value A(jito), we 
combine the two RGEs of Eq. (lilt to obtain 



A 6 (ju) =X h {lio) 



c(fl. s (Mo)) ' 



(13) 



with 



c{a) — a c ° ^ 1 + (cj — b\Co) a 
1 



[(ci - b\c a ) 2 + c 2 - b\c\ + b\c Q - b 2 co] a 2 



c . = t_ h .= h 
C, -]3o' '"ft' 



(14) 



2 Note that only the e-print of Ref. l6lll discusses the renormalization of the 
pseudoscalar cun'ent. 



3 The numerical value of a s (Mz) has to be set in accordance with the parton 
sets that are used, see below. 
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FIG. 5: The cross section a(pp — » (bb)H + X) (in picobams) at 
(a) LO, (b) NLO, (c) NNLO for the LHC. The axes labels are F = 
\og w {n F /M H ) and R = \og l0 (n R /M H ). Thus, the point )X R = M H , 
[if = 0.25M// corresponds to R = 0, F = —0.6. The Higgs boson 
mass is set to Mh = 120 GeV. 
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FIG. 6: The cross section cs{pp — > (bb)H +X) (in picobarns) at (a) 
LO, (b) NLO, (c) NNLO for = 1 .96 GeV. The notation is the same 
as in Fig.|5| The Higgs boson mass is set to Mh =120 GeV. 



suits in the physical hadronic cross section. The variation of 
the hadronic cross section with and jif> is thus an indication 
of the size of higher order effects. 



j3™ and y* have been defined in Eq. (TT21 . Both a s (fi) and 
a s (/Uo) are calculated from a s (Mz) using the procedure de- 
scribed above. Working at LO (NLO, NNLO), we truncate the 
term in braces at order a (a 1 , a 2 ). 

Convolution of the partonic cross section with the parton den- 
sities cancels the jj, F dependence up to higher orders and re- 



IV. RESULTS 

The analytic expressions for the partonic cross section are 
quite voluminous and will be deferred to the appendix. In 
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FIG. 7: Cross section for pp -> (bb)H+X at V* = 14TeV, M H = 
120 GeV. (a) /x« dependence for flf = 0.25M#; (b) flf dependence 
for fi R =M H . 
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FIG. 8: Cross section for pp -> (bb)H+X at y/s = 14TeV, M H = 
300 GeV. (a) fl R dependence for flf = 0.25M#; (b) flf dependence 
for Hr = M h . 



this section, we study the behavior of the NNLO result with 
respect to variations of the input parameters, in particular the 
Higgs boson mass and the collider type (LHC and Tevatron). 
Special emphasis is placed on the variation of the results with 
the renormalization and factorization scale, from which we es- 
timate the theoretical uncertainty of the prediction for Higgs 
boson production in bb annihilation. 

Because the cross sections for the neutral Higgs bosons in 
bb annihilation differ only in the magnitudes of the Yukawa 
couplings (within our approximations), we will restrict our 
discussion to the production of a Standard Model Higgs bo- 
son. In the limit that supersymmetric partners are heavy, their 
virtual contributions are insignificant and the predictions for 
supersymmetric Higgs bosons can be obtained from the Stan- 
dard Model values by rescaling them with the proper coupling 
constants (cf. Eq. (I10» . 



All the numerical results have been obtained using 
Martin-Roberts-Stirling-Thorne (MRST) parton distribu- 
tions. In particular, we use the MRST2001 sets at 
LO [a s (M z ) = 0.130] and NLO [a s (M z ) = 0.119], and 
MRSTNNLO jH at NNLO \a s (M z ) =0.1155]. 

In order to obtain an overall picture of the renormalization and 
factorization scale dependence of the cross section, we plot 
a(pp — > (bb)H +X) for = 14TeV as a function of the 
two parameters fip and jj.r in Fig.|5] The corresponding plot 
for the Tevatron, i.e., (j(pp -> (bb)H+X) for y/s = 1 .96 TeV, 
is shown in Fig.[6] The Higgs boson mass is fixed to Mh = 
120 GeV. Subpanels (a), (b), and (c) show the LO, NLO, and 
NNLO prediction, respectively. Note the extremely large vari- 
ation of the scales, by a factor of 10 above and below Mh- 
Apart from the region of large flf and small fig, one observes a 
clear reduction of the scale dependence with increasing order 
of perturbation theory, both for flf and fiR. Notably, we find 
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FIG. 9: Cross section for pp^> (bb)H + X at ^fs = 1 .96 TeV, M H = 
120 GeV. (a) fl R dependence for flf = 0.25M#; (b) flf dependence 
for fi R =M H . 



minimal radiative corrections and a particularly weak depen- 
dence on the renormalization and factorization scales in the 
vicinity of (fi R ,^i F ) = (Mq,0.25M$) = (p.R,p.f). This agrees 
with the observation of Ref. lfl7ll that the proper factorization 
scale for this process should be around jU/r = /A. 

To illustrate this observation, we display separately the /Ir- 
and (if -variation of the cross section at the LHC in Fig.Qfor 
M H = 120 GeV, and in Fig. HO for M H = 300 GeV. In subpan- 
els (a), the factorization scale is fixed to fip = p,p = Q.25Mh, 
and the renormalization scale is varied within 0. 1 < PLr/Mh < 
10. In subpanels (b), the renormalization scale is fixed to 
jUr = P-r — Mh, and the factorization scale is varied within 
0.1 < [If/Mm < 10. The reduction in the scale dependence 
with increasing order of perturbation theory is clearly visible. 
As opposed to the lower order curves which have a monotonic 
dependence on Hg/p within the displayed range, the NNLO 
curves develop a maximum, so that it is possible to define a 
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FIG. 10: Cross section for Higgs boson production in bottom quark 
annihilation at (a) the LHC and (b) the Tevatron (Run II) at LO 
(dotted), NLO (dashed) and NNLO (solid). The upper (lower) line 
corresponds to a choice of the factorization scale of flf = 0.7M# 
([if — 0.1 Mu). The renormalization scale is set to fi R = Mfj. 



"point of least sensitivity" for them. In all cases, this falls 
nicely into a region where the radiative corrections are small. 
Note also that the central values for the NNLO curves are per- 
fectly consistent between panels (a) and (b). These observa- 
tions confirm that jXp = 0.25 Mh and /1r = Mh are indeed the 
appropriate scale choices for this process. 

The corresponding curves for Run II at the Tevatron are shown 
in Fig.|9] (we only show results at the Tevatron for Mh = 
120 GeV). As opposed to the LHC, the reduction of the renor- 
malization scale dependence with increasing order of pertur- 
bation theory is less drastic. One even observes a slight in- 
crease in the /Ir dependence between NLO and NNLO. How- 
ever, the absolute variation is very small. The factorization 
scale dependence is quite similar to that observed for the LHC. 
Again, the central values for the NNLO curves of subpanels (a) 
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and (b) coincide nicely. Note that the cross section at the Teva- 
tron is typically about two orders of magnitude smaller than 
at the LHC. 

Fig. llOf a) shows the LO, NLO, and NNLO predictions for the 
cross section o(pp — > (bb)H + X) at the LHC as a function 
of the Higgs boson mass. The two curves at each order cor- 
respond to two different choices of the factorization scale, 
ji F = 0AM H and Hf = 0.7 M H - From subpanels (b) of Figs.0 
and[8]one can see that this roughly defines the maximal pif- 
variation at NNLO between 0.1 M H and 10M H - Since the 
renormalization scale dependence is very weak [cf. subpanels 
(a) of Figs.0and[8], we fix Hr = £ir = Mh- Taking the width 
of these bands as an indication of the theoretical uncertainty, 
we observe an improvement of the accuracy of the prediction 
for M H = 120GeV from 70% at LO to 40% at NLO to 15% 
at NNLO. At larger Higgs boson masses, the scale uncertainty 
is smaller, amounting to 40% at LO, 17% at NLO, and 5% at 
NNLO for M H = 300 GeV. 

The cross section for the Tevatron at y/s = 1 .96 Te V center- 
of-mass energy is shown in Fig. llOr b). Here the renormaliza- 
tion scale dependence within the range 0.1 < JJ-r/Mh < 10 
at NNLO is larger than the factorization scale dependence 
(cf. Fig.|9j- Nevertheless, we apply the same prescription as 
for the LHC and plot the LO, NLO, and NNLO cross section 
at (jUs,Mf) = (l,0.1)Afe and {}1r,IIf) = (1,0.7)M H . This 
is justified since /i^-variation on absolute scales is still very 
small, in particular if it is restricted to a more reasonable 
range of about a factor of five above and below Mh- We ob- 
tain a reduction of the scale uncertainty at the Tevatron for 
M H = 120 GeV from around 60% at LO, to 30% at NLO, to 
10% at NNLO. 



V. CONCLUSIONS 



We have computed the total cross section for Higgs boson 
production in bb fusion at NNLO in QCD. We have argued 
that the NNLO plays an exceptional role in this process, as 
it incorporates all subleading logarithms at order a s 2 . The 
results are very stable with respect to changes of the renor- 
malization and factorization scales. We find that the radia- 
tive corrections are particularly small at factorization scales 
of around Uf = M^jA, in agreement with the arguments of 
Refs. E3I21. 

We conclude that the inclusive cross section for Higgs boson 
production in bottom quark annihilation at hadron colliders is 
under good theoretical control. 
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APPENDIX A: PARTONIC RESULTS 



It is convenient to write the partonic cross section in the fol- 
lowing way: 



&ij (x) = (7° Ay 0) , i,j e{b,b,g,q,q}, 



(Al) 



where Ojj is the cross section for the process ij —> (j>+X. i and 
j label the partons in the initial state, (j) means either a scalar 
or pseudoscalar Higgs boson, and X denotes any number of 
quarks or gluons in the final state. Here and in what follows, 
q denotes any of the light quarks u,d,s,c. The normalization 
factor, Co, is 



12 Ml 



(A2) 



The correction terms are written as a perturbative expansion: 

A l7 (*) = A|°) (x) + £ A« ft + (f ) 2 A[f(x) + 6{a]) . 

(A3) 

Explicit dependence on the number of active flavors n/ ap- 
pears only at NNLO. Because the terms are large and cumber- 
some, it is convenient to write 



(A4) 



All results will be presented for the scale choices /J.f = /J.r = 
Ma . The corresponding expressions for general values of jj,f 
and Hr can be reconstructed from renormalization scale in- 
variance of the partonic, and factorization invariance of the 
hadronic cross section. 4 



The analytic results including all scale dependences can also be obtained 
from the authors upon request. 
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1. The bb subprocess 

In the VFS approach, the tree-level bb annihilation term is the LO contribution. Thus, this is the only term for which (x) does 
not vanish. The LO contribution to bb — > <j> +X is 

A$(x) = S(l-x). (A5) 

The NLO contribution is 

4-8C2C-/, x 16 „ „ \ 16 + 8x + 8x 2 , Ax-Ax 2 8 In fx) 8+4x + 4x 2 , . 
A S W= ^5(l-x) + T ^(l-x) ln(l-x) + 3 3 + 3 ln « > (A6) 



where @ n (l—x) 



ln"(l-jc) 



, and C2 = ^ 2 /6 ~ 1.64493, £ 3 ~ 1.20206. 



1-X 

At NNLO, the contributions are (£ 4 = 7r 4 /90 w 1 .08232) 

(2)a 115 + 116^-156^-19^ 404-396^-1146^ 

lb = is 8(1 -x) 9o(l-x) 

204- 200 C 2 , ^ 44 128 
+ r ^$> l (l-x)--$> 2 (l-x) + — &{l-x) 

l28 + 64x + 64x 2 3/l , 140 + 40x + 92x 2 - 8x 3 , 2 
g In (1 -x) + g In 2 (1 -x) 

248 In 2 (1-x) In (x) 248 + 168x + 168x 2 , , „ 
— + In (l-x)ln(x) 

604+138x + 423x 2 + 44x 3 . 200+ 100x + 100x 2 „ , 
27 In (1 -x) + C2 In (1 - x) 

„ lnfl-x)ln(x) 216+110x+160x 2 -24x 3 , , . 

+ 2A-— — - — — ln(l -x) ln(x) 

1 — X y 

148 ln(l -x)ln 2 (x) 148 + 110x+ 110x 2 , N ~ 

+ — - K -L In 1 - x In 2 x 

9 1 —x 9 

201n(l-x)Li 2 (l-x) 20-78x-78x 2 . T . 

+ T fcP 9 ln(l-x)Li 2 (l-x) 

4640+1017x + 2958x 2 + 721x 3 140 - 15x+ 147 x 2 - 8x 3 „ 
+ 3^ ~ 9 C2 (AT) 

382 + 191x + 191x 2 146 In (x) 164 £ 2 In fx) 23 In 2 (x) 44 In 3 (x) 

9 ^ ~ ~9~ 1 -x + "9" 1 -x ~ T 1-x ~ 27 1 -x 

876 + 249x + 444x 2 + 38x 3 , , , 164+ 126x+ 126x 2 „ , , . 
+ 54 In (x) C2 In (x) 

138+73x+ 115x 2 - 12x 3 , , 44 + 21x + 21x 2 + 4x 3 , , . , 

H In fx) H In (x) 

18 w 27 y ' 

Li 2 (l-x) 58 Li 2 (l -x)ln(x) 142Li 3 (l-x) 64 . ( 1-x 
+ 1-x + ~9 ~ x 9 l - x 9 (1-x) Ll3 

7-51x-10x 2 +10x 3 T . ^ . 58 + 88x + 88x 2 + 2x 3 T . „ . . . 
H ^ Li 2 ( 1 - x) Li 2 ( 1 - x) In (x) 

- X - lx2 - 2x ' Li 2 (1 -x 2 ) + j Li 2 (1 -x 2 ) In (x) 

142 + 37x + 37x 2 + 6x 3 . , 64 + 94x + 94x 2 -6x 3 T . ( 1-x 
H Li 3 (1 -x) H Li 



9 " JV " "' 9 " 3 V x 



— Li 3 (l-x)--Li 3 



18 " J v " " ' 18 " J V x 2 J 3 



/l-x\ / 1-x 

L13 h-r- - Ll 3 - 



\l+x/ V 1+* 



12 



.(2)1? 2-10&+18& , , 56-72& a /i x , 8 » n 

8+4x + 4x 2 , , . 40 + 8x + 32x 2 , . 16 lnfl -x) \n(x) 
n ln 2 (l-x) + ln(l-*)- 



- - (A8) 

16 + 8x + 8x 2 . , lOlnfx) 21n 2 (x) 2Li 2 (l-x) 56+x + 55x 2 

lnfl -jc) lnfx) H — H — — 

9 v ' w 9 1-x 3 1-x 9 1-x 81 

8 + 4x + 4x 2 10 + 3x + 7x 2 . . 12 + 7x + 7x 2 , 2 . . 2_... , 
C2 In (x) — In 1 (x) + - Li 2 (1 -x) . 



9 9 w 18 

Note that the @ n terms in this result could also be derived by other methods Ol64ll65ll65l67ll68ll69ll . 



2. The bg subprocess 



The bg <j> +X subprocess first enters at NLO, where the contribution is 



At NNLO, the contribution is 

A2)A _ .(2)A _ 
^bg ~ % - 

257x-514x 2 + 514x 3 



x — 2x + 2x 



In (1-x) 



3x-10x 2 + 7x 3 x-2x 2 + 2x 3 



lnfx) 



+ ■ 



144 

16-59x + 272x 2 -237x 3 
16 

16 + 28x-731x 2 + 726x 3 



In 3 (1-x) 

2/1 



llx-94x 2 + 62x 3 , , 



48 

65x-508x 2 + 774x 3 



ln(l-x) 



ln(l-x)lnfx)- 



35x 


8 

-70x 2 - 


h70x 3 




12 




3x- 


1- 174x 2 


-98x 3 , 



24 

3x + 6x 2 + 6x 3 



In 2 (1-x) In (x) 



In (1-x) ln z (x) 



24 

77x+ 134x 2 — 86x 3 _i.s — u.i - r u.\ 

— ln(l -x) Li 2 (1 - x) + ln(l - x) Li 2 (l -x 2 ) 

208-411x- 1350x 2 + 1781x 3 16 - 31x+ 176x 2 - 169x 3 „ 161x- 322x 2 + 322x 3 



288 

32x + 536x 2 -993x 3 
48 

47x-604x 2 +1028x 3 
96 



16 



48 



, , . 7x-230x 2 +134x 3 - , , , 
In W + 24 & W 

, 2 , . 35x+ 146x 2 +12x 3 , 3 . , 

In fx) In fx) 

w 144 W 



48-43x+152x 2 + 236x 3 T . . 10x + 34x 2 - 17x 3 r . „ . . 

H — Li 2 (1 — x) Li 2 (1 — x) ln(x) 



24 

19x + 40x 2 + 21x 3 



48 

14x + 8x 2 + 9x 3 



, . 2 , 9x+18x 2 + 10x 3 ~ 
Li 2 (1 -x ) — Li 2 (1 -x ) lnfx) 



12 



r . , 13x+118x 2 -18x 3 
L13 ( 1 — x) H 77; L13 



9x+18x 2 -14x 3 
48 

3x + 6x 2 + 6x 3 r 



Li 3 (l-x 2 ) + 



12 

9x+18x 2 -14x 3 



U.(£|-L„|- 



48 

1 -x 
1 +x 



-Li 3 



1-x 



1 -x< 



(A9) 



(A10) 
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3. The bq subprocess 

All remaining components contribute only at NNLO and beyond. The contribution to bq —* <j) +X, where q is a light (u,d,s or c) 
quark, is: 

A2)A J2)A J2)A ,(2)A 4 + 3x - 3x 2 - 4x 3 2 

< = A L/ = A W =4, = 9 ln 2 (l-x) 

+ 2x + 2x 2 ln2(1 _ x) ln(x) + 4 - 57x +g x2 - 22x3 ln(1 _ x) 

x _4x 2 -4x 3 2x + 2x 2 
ln(l -x) \n(x) ln(l -xj In (x) 

4x + 4x 2 „ ST . ^ . 208-915x+1410x 2 -703x 3 
+ 3 1" (1 Li 2 (l-x) — 

4 + 3x-3x 2 -4x 3 „ 93x-264x 2 + 20x 3 , , . 2x + 2x 2 „ , , . (A11) 

9 ^ + Tos ln (x) 3— Cz ln {x) 

3x+15x 2 + 40x 3 , . x+x 2 ,,,, 16-3x + 21x 2 + 8x 3 r . ^ 
72 ln to + "^6~ ln to + ig Ll 2 (! -*) 

- (x+x 2 ) Li 2 (1 - x) ln (x) - 2 JL±1?L Li3 (1 -x) + Li 3 {^~) > 

t$ F = 6& F = =0. 

»</ fcg foj 



4. The gg subprocess 

The contribution to gg — ► +X is: 

Ag )A = - (x + 2x 2 - 3x 3 ) In 2 (1 - jc) - X + 4 ^ + 4x ' In 2 (1 -x) ln (x) 

23x + 52x 2 -75x 3 , , 5x+16x 2 -4x 3 , . . . 

H ln(l -x)H ln(l -x) ln(x) 

8 4 



x + 4x 2 + 4 



3 



+ ~ ' ' " X In ( 1 - x) In 2 (x) - (x + 4x 2 + 4x 3 ) ln ( 1 - x) Li 2 ( 1 - x) 

163x+ 1528x 2 - 1691x 3 - „ , „ „ 
j2g + (* + 2^ 2 - 3x 3 ) C 2 

54x + 312x 2 -223x 3 , , . x + 4x 2 + 4x 3 „ , , . 
g4 1" to + ^ C2 In (x) 

16x+lllx 2 -43x 3 , 2 , , 7x + 25x 2 + 34x 3 , N 
64 ln « + 48 ln « 

- 4 - 15 - 2 ' 62x3 u 2 (i - x) + llx+4 t 2 + 3 ° x3 ii 2 (i - x) m to (M2) 

16 16 

x 2 -6x 3 . 2 , 3x + 6x 2 + 38x 3 , 2 , 
+ — ^ — Li 2 (l-x )+ — Li 2 (l-x )ln(x) 

x + 3x 2 + 18x 3 . „ . 15x + 60x 2 + 30x 3 T . / l-x 
+ g Ll 3 (1 -x) — L13 



5x+10x 2 + 74x 3 r . - 7N 3x + 6x 2 + 70x 3 T . / l-x 2 
-Li 3 (l-x) — L13 3- 



128 v > 128 



x + 2x 2 + 2x 3 
32 



L,(^)-Li3 



( 1+*) 



Ag )f = 0. 
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5. The bb subprocess 

The contribution to bb — > <j) +X is: 



A<f = A[f = 8 + fa -f- 8 - 3 In 2 (1 - + 1^ In 2 (1 - ,) 1„ (,) 

8- 138x+174x 2 -44x 3 , , lOx - 20x 2 - 24x 3 , „ . . 

+ 27 ln(l-x) ln(l-x) ln(x) 

4 ln(l-x)ln 2 (x) 16 ln(l -x) Li 2 (1 - x) 8 ln(l -x) Li 2 (l -x 2 ) 

+ 9 f+JC ~9~ T+x + 9 f+~x 

wi \i 2/ \ j 16 + 16x + 32x 2 
ln(l -x) hr (x) H ln(l -x) Li 2 (1 -x) 



4+10x+14x 2 ^ N , 16 + 16x + 32x 2 
In ( 1 — x 1 lr- 

;-4x + 4x 2 , „ w . 2 s 52-357x + 510x 2 -205x 3 



ln(l-x)Li 2 (l-x 2 )- 



9 ' v y 81 

8+6x-6x 2 -8x 3 „ 117x-279x 2 + 20x 3 , , N 4x + 4x 2 „ , , N 
^ C2 + ^ In (x) £> ln(x) 

x+llx 2 + 34x 3 , 2 , . 28- 12x+ 17x 2 , . 14 In 3 (x) 

36 ln « + 54 ln W - 27 T+T 

4 Li 2 (l-x)ln(x) 8 Li 2 (l-x 2 )ln(x) 4Li 3 (l-x) 4 / l-x\ ( A1 3) 

+ o r~; n o — r~ 1~ ttt~ — 7L13 -- 



3 1+x 9 1+x 3 1+x 3(l+x) V x 



lLi 3 (l-x 2 ) 1 / l-x 2 \ 20 

+ ^ 7~ - + »/,, n Li 3' 



3 1+x 9 (1+x) V * / 9 (1+x) 



, 1 — x \ ( 1 — x 
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16-12x+16x 2 +llx 3 r . 12 + 8x + 29x 2 y . 
H jj Li 2 (l -x) Li 2 (l -x)ln(x) 

x T . u 2N 16-llx+llx 2 T . , . . 12-27x + 5x 2 T . 
+ -Li 2 (1 -x 2 ) + Li 2 (1 -x 2 ) ln(x) + L13 (1 -x) 

12- 16x-9x 2 -2x 3 / l-x\ 12-13x+13x 2 , 2 n 
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/ 1 — x\ . ( 1 —x 



9 J V x J 36 
x + 3x 2 T . / l-x 2 \ 20- 13x+13x 2 L . 
36 Ll3 {- — ) 9 Ll3 



1 +x 



6. The subprocess 

The contribution to qq — ► +X is e {m,c/,.s,c}): 



^ : 2 :. .4 2x-8x 2 + 6x 3 x-2x 2 -3x 3 

v-3 1 „3 



<r= - - Q — m*)+-5-^w 

-" 9 — ^v- -v ■ 9 



4x J 2x J 

— Li 2 (l-x) + — Li 2 (l-x 2 ) , (A14) 



4f = °- 
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